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system. Making use of the coupling between the heat bath and the external modulation 
as a small perturbation we construct a Langevin equation with multiplicative noise 
and space dependent dissipation and the corresponding Fokker-Planck-Smoluchowski 
equation in the overdamped limit. We examine the thermodynamic consistency 
condition and explore the possibility of observing a phase induced current as a 
consequence of state dependent diffusion and, necessarily, nonlinear driving of the 
heat bath by the external noise. 
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In recent times the phenomena of noise induced transport under nonequihbrium 
condition have gained wide interdisciphnary interest where the interplay of fluctuations 
and nonlinearity of the system plays an important role [H El El HI [5] . Exploitation of 
the nonequihbrium fluctuations present in the medium helps to generate phase induced 
directed motion of the Brownian particle. Presence of spatial anisotropy in the potential 
together with nonequihbrium perturbations enables one to extract the useful work from 
random fluctuations without violating the second law of thermodynamics [I3j. This 
lead to its wide applicability in explaining the dynamics of molecular motors [21 [6], 
directed transport in photovoltaic and photoreflective materials [7] , and the efficiency of 
tiny molecular machine in a highly stochastic environment [3, HI [8] , realization of ratchet 
effect in cold atom [9], construction of artificial molecular rotors that produce controlled 
directional motion mimicking molecular motor protein [10]. In some special cases one 
can generate directed motion even in a symmetric potential due to state dependent 
diffusion. For such systems the state dependent diffusion coefficient, -D(g), felt by the 
Brownian particle could arise either due to space dependent friction or presence of local 

hotspots [21 El nulls]. 

To the best of our knowledge, in almost all the above mentioned cases the 
corresponding Langevin equation is either written phenomenologically or is constructed 
from a microscopic system heat bath Hamiltonian model where the associated heat bath 
is in thermal equilibrium. To generate directed motion one then applies an external 
perturbation, time periodic force or correlated random force, to break the symmetry 
of the force field as equilibrium thermal fiuctuations are unable to create spontaneous 
symmetry breaking. Symmetry breaking can be also achieved by considering a nonlinear 
coupling between the Brownian particle and thermal heat bath thereby generating 
a multiplicative noise term in the resulting Langevin equation which in turn gives 
a position dependent diffusion term effectively creating a phase induced bias in the 
dynamics. In the present paper, we propose a system heat bath model where the 
heat bath is weakly modulated by an external noise. Although the microscopic model 
we present here has a close kinship to our earlier approach [I3] to study escape 
from a metastable state within the context of external noise modulated heat bath, 
the present formalism differs from our earlier model in the following way. The heat 
bath-external noise coupling is considered to be nonlinear and in addition to that the 
system is also nonlinearly coupled with the heat bath thereby resulting in a nonlinear 
multiplicative Langevin equation and the corresponding Fokker-Planck-Smoluchowski 
equation with space dependent diffusion. We then explore the possibility of observing 
directed transport as a result of phase difference between the coupling function and the 
periodic potential in which the Brownian particle is moving. Our theoretical model can 
be tested experimentally to study the directional motion of artificial chemical rotors in 
photoactive solvent [lO] . To observe the effect of external stochastic modulation one can 
carry out the experiment in a photochemically active solvent (the heat bath) where the 
solvent is under the influence of external monochromatic light with fluctuating intensity 
of a wavelength which is absorbed solely by the solvent molecules. As a result of it the 
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modulated solvent heats up due to the conversion of light energy into heat energy by 
radiationless relaxation process and an effective temperature like quantity develops due 
to the constant input of energy. Since the fluctuations in the light intensity result in the 
polarization of the solvent molecules, the effective reaction field around the reactants 
gets modified [H]. 

To start with we consider a classical particle of unit mass bilinearly coupled to 
a heat bath consisting of mass weighted harmonic oscillators characterized by the 
frequency set ujj. In addition to that the heat bath is nonlinearly driven by an external 
noise e(t). The Hamiltonian for the composite system is 

H = Hs + Hb + HsB + Hint-, (1) 

where Hg = + V{q), is the system's Hamiltonian with q and p being the 

coordinate and momentum of the system particle, respectively, and V{q) is the potential 
energy function of the system. Hb + Hsb = Ef=i[{p'j/'^) + KV2){a;j - Cj/(g)}2] 
where {xj,pj} are the variables for the j-th bath oscillator. The system-heat bath 
interaction is given by the coupling term CjUJjf{q) with cj being the coupling strength. 
We consider the interaction. Hint = J2f=i f^jg{xj)^{'t), between the heat bath and 
the external noise e{t) where Kj denotes the strength of the interaction and g{xj) 
is an arbitrary analytic function of the bath variables, in general nonlinear. This 
type of interaction makes the bath variables explicitly time dependent. A large 
class of phenomenologically modeled stochastic differential equation may be obtained 
from a microscopic Hamiltonian for particular choice of coupling function g{xj). In 
what follows we have chosen g{xj) = x'^/2, which makes the spring constants of the 
bath oscillators time-dependent. The external noise is stationary, Gaussian with the 
properties (e(t))e = and (e(t)e(t'))e = 2D5{t — t'), where D is the strength of the 
external noise and (■)e implies averaging over the external noise processes. From Eq.([l]) 
we have the dynamical equations for the system and bath variable 

q(t) = -V'{q{t)) + f{q{t)) ^ c,u;]{x,{t) - c,/(g(t))}, (2) 

j 

x,{t) + {uj] + ^,€{t)}x,{t) = c,u^f{q{t)), (3) 

where we have used g{xj) = x^/2. To solve Eq.(l3]) for Xj we assume a solution of the 
form 

Xj{t) = x°j{t) + i^jx]{t), (4) 

where x^(t) is the solution of the unperturbed equation of motion 

x%t)+u;]x%t) = c,u;^f{q{t)). (5) 

The physical situation that has been addressed here is the following, we consider that 
at t = 0, the heat bath is in thermal equilibrium in the absence of the external noise 
e{t). At t = 0_|_ the external noise agency is switched on and the heat bath is modulated 
by e{t) [I3]. Then x^(t) must satisfy the equation 

x]{t)+uj]x]{t) = -x'^{t)e{t), (6) 
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Uj Jo 



dt' smujj{t - t')Xj{t')e{t'). 



(7) 



The formal solution of Eq. ([5]) is given by 



= a;°(0) cosujt + sincj^t + CjUJj / dt' smujj{t - t')f{q{t')) (8) 



where x°(0) and Pj'(O) are the initial position and momentum, respectively, of the jth 
oscillator. Now using this solution in Eq.Q we have, after an integration by parts, the 
equation of motion for (t) which gives the equations of motion for the bath variables 
Xj{t) (from Eq.(jl])) as 



x,{t)= xm-c,f{qm 



COS ujd H — 



UJj 



sin ujjt 



UJj Jo 



dt' COS ujj{t-t')f'{q{t'))q{t') 
dt' sinujit - t')e{t')x%t'). 



(9) 



Using the above solution in Eq.([2]) we finally obtain the equation of motion for system 
variable as 



m = - V'{q{t)) 



|x^(0) — Cj/(g(0))| cosujjt H — ^ sinujjt 

- E ll dt' cos^,(t - t')f'{q{t'))p{t') 

r-t 

~ f'{Q{t)) E ^i'^j^j J sinco'j(t — t')e{t')x^j{t'), 



(10) 



where p{t) = q{t) is the generalized momentum of the system variable. This equation 
can be rewritten as 



q{t)= -V'{q{t))-f'{q{t)) dt'^{t - t') f'[q{t')]p{t') + f {q{t))F{t) 



- f\<lit)) E (^j^j^j / dt' sinujjit - t')x]{t')e{t'), 



(11) 



where we have defined 7(t) and F{t) as, 7(t) = Yl!j=iC^j^'^ cosujit) and F{t) = 
J2f=iCjUj'j[{xj{0) — Cj f{q (0)} cos ujjt + {pj{0) /cuj) sinujjt]. At this point, we note that, 
the forcing term F{t) is deterministic as expected. It ceases to be deterministic when 
it is not possible to specify all the x°(0)-s and p°(0)-s, i.e., the initial conditions of all 
the bath variables, exactly. The standard procedure to overcome this problem is to 
consider a distribution of a;°(0) and p°(0) to specify the statistical properties of the bath 
dependent forcing term F(t). The distribution of the bath oscillators is assumed to be 
a canonical distribution of Gaussian form 

Hb + Hsb' 



W 



a;°(0),p°(0) =Z-'exp 



(12) 
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where Z is the bath partition function. This choice of the distribution function of bath 
variables makes the internal noise F{t) Gaussian. It is now easy to verify the statistical 
properties of F{t) as {F{t)) = and {F{t)F{t')) = 2kBT-f{t - t'). where ks is the 
Boltzmann constant and T is the equilibrium temperature. (■) implies the average over 
the initial distributions of bath variables which is assumed to be a canonical distribution 
of Gaussian form as given in Eq. f|T2|) . The second relation is the celebrated fluctuation- 
dissipation relation [15] which ensures that the bath was in thermal equilibrium at t = 0. 

To identify Eq.l ITT]) as a generalized Langevin equation we must impose some 
conditions on the coupling coefficients Cj and kj, on the bath frequencies Uj and on 
the number of the bath oscillators that will ensure that 7(t) is indeed dissipative 
and the last term in Eq. ffTTl) is finite for oo. A sufficient condition for 7(t) 

to be dissipative is that it is positive definite and decreases monotonically with time. 
These conditions are achieved if A^ — >• cx) and if cju'j and Uj are sufficiently smooth 
function of j [17]. As A^ — oo one replaces the sum by an integral over u weighted 
by a density of state p{uj). Thus to obtain a finite result in the continuum limit the 
coupling function q = c{u) and = are chosen [131 [IB] as c{ij) = c^juj^ ^ and 
K[io) = kq, where cq and kq are constants and Tc is the correlation time of the heat 
bath. The choice k,{uj) = kq is the simplest one where we assume that every bath 
mode is excited with the same intensity. This simple choice makes the relevant term 
finite for N ^ oo. Consequently ■y{t) becomes, 7(t) = (cI/tc) J du p{uj) cos out. 1/tc 
may be characterized as the cutoff frequency of the bath oscillators. The density of 
modes p^uj) of the heat bath is assumed to be Lorentzian, p{uj) = (2/7r)[rc/(l + uj^t^)]. 
The above assumption resembles broadly the behavior of the hydrodynamical modes 
in a macroscopic system [18]. With these forms of p{uj), c{uj) and k{uj) we have the 
expression for 7(t) as 7(t) = (cq/tc) exp(— t/rc), which reduces to 7(t) = 2cl6{t) for 
vanishingly small correlation time and consequently one obtains a 5-correlated noise 
process. 

Taking into consideration of all the above assumptions and assuming that the 
system variable evolves much more slowly in comparison to the external noise e{t), 
in the limit 0, Eq. fllll) reduces to 

q{t) = -V'{q{t)) - l[nq{t))fm + /'(^W)^^^ + l^of{q{t))nq{t)W), (13) 

where 7 = Cq is the dissipation constant and the Langevin force F{t) is characterized 
by the statistical properties 

{Fit)) = 0, {F{t)F{t')) = 2^kBT5{t - t'). (14) 

Eq.f ll3l) is the generalized Langevin equation for the system variable. At this juncture, 
it is noteworthy that for /(g) = g, Eq.(IT3ll reduces to q{t) = —V'{q(t)) — -yqit) + F(t) + 
7/€og(t)e(t). Thus for linear system-bath coupling (i.e. for f{q) = q), our Hamiltonian 
given by Eq.([T]) may be the starting point for the construction of a Langevin equation 
with both additive and multiplicative noise, which has numerous applications in various 
field of physics, e.g. phase transition etc. For harmonic potential, this equation has 
been extensively studied by many authors in various contexts |19] . 
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In the above Langevin equation fll3p the noise terms (internal and external) appear 
multiplicatively and the dissipation is space dependent. Using the method of van 
Kampen [19j for nonlinear stochastic differential equations the Fokker-Planck equation 
corresponding to the Langevin equation (fT3|) is given by [131 [IH] 



+ {\(q)kBT + ',\ID 0^ (15) 

where P = P{q,p, t) is the phase space probability density function and \{q) = 'y[f'{q)Y 
is the space dependent dissipation function. Instead of handling two noise processes 
(internal and external) independently, one can define an effective noise process ^{t) and 
an auxiliary function G{q) to obtain the same Fokker-Planck equation f|T5l) from the 
following Langevin equation 

q = -V'iq)~\iq)q + G{q)^{t), (16) 

with 

m) = o,mm) = mt-t'), (ir) 

m = llf'iq)?, G{q) = f\q)^JlkBT + D{^K,YP{q). (18) 

That the Langevin equation ( |T6i) gives the same Fokker-Planck equation ( |T5l) can 
be verified by using van Kampen's methodology [19]. The construction of Langevin 
equation using an effective noise term and an auxiliary function has been done earlier 
in the configuration space by Wu et al |20], whereas we have written the Langevin 
equation ( fT6|) in the phase space. Thus as far as the equation for the evolution of 
probability density function is concerned, Eq.(|T6l) is the equivalent description of the 
stochastic differential equation (fT3l) . Eq. (fT6|) is one of the key results of this work as it 
incorporates the effects of thermal noise F{t) and the external noise e{t) in an unified 
way even when the underlying noise processes are multiplicative due to the nonlinear 
system-bath coupling and nonlinear modulation of the heat bath by an external noise. It 
is important to mention here that Eq. (fT6!l describes a thermodynamically open system 
where there is no fluctuation-dissipation relation so that the system will not reach at 
usual thermal equilibrium, instead, a steady state is attainable for large t [I3l [19] . From 
the computational point of view, generation of a single multiplicative noise process is 
much more economical than to generate two separate multiplicative noise processes. 

In Eq.( !T6|) the noise is multiplicative and the dissipation is space dependent. In the 
case of large dissipation one eliminates the fast variables adiabatically to get a simpler 
description of the system dynamics. The traditional approach to the elimination of fast 
variables for multiplicative noise processes does not always give the correct description. 
In order to get the correct Langevin equation in the overdamped limit we follow the 
method of Sancho et al ^21j and then using van Kampen's lemma [2^ and Novikov's 
theorem [23] we get the Fokker-Planck-Smoluchowski equation corresponding to Eq. (|T6il 
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for the probability density P{q,t) in the configuration space 



dP{q,t) d 1 



dt 



dqX{q) 



V'iq) + 



d G\q) 



P{q,t). 



(19) 



dq X{q) 

In the ordinary Stratonovich description [21] the Langevin equation corresponding to 
the Fokker-Planck equation ( fTOll is 

V'{q) G{q)G'{q) , G{q) 



(20) 



The above equation ( l20l) differs from the Langevin equation, obtained by using the 
traditional way of adiabatic elimination of fast variable, due to the presence of the 
second term on the right hand side. This term, G(g)G'(g)/[A(g)]^, represents the effect 
of multiplicative noise in the process of elimination of fast variable |21j . 

The stationary solution of Eq. fllQI) contains inhomogeneous effective temperature 
like term, a generalization of Boltzmann factor for state dependent diffusion in 
open system, which arises due to the entanglement of the external driving with the 
nonlinearity of the system heat bath coupling, a well known effect in several contexts, 
e.g., Landauer Blow torch effect In the absence of external bath modulation, 

i.e., when G{q) 
function, Peg{q) 

the overdamped limit we then have the stationary current as 



- f'{q)\^^ykBT, Eq. (fT9ll gives the correct equilibrium distribution 
JV exp[—V{q)/kBT] with A/" being the normalization constant. In 



J 



A(g) 



V'iq) + 



d (G\q)- 
dq I A(g) 



(21) 



Integrating the above equation we have the expression of stationary probability 
distribution in terms of stationary current 



where 



PsM) 



-4>{q) 



G\q)/\{q) 
V'{q') 



G\q')l\{q') 



G'^(O) 
A(0) 



dq' 



Pst{0)-J r X{q')e'''^'''^dq' 
Jo 



V'iq') 



ksT + D^KlPiq' 



-dq'. 



(22) 



(23) 



We then consider a symmetric periodic potential with periodicity 27r, V{q) = V{q + 27T) 
and the periodic coupling function with the same periodicity as the potential, /(g) = 
/(g + 27r). Now applying the periodic boundary condition on Pst{q), Pst{q) = -Pst(g + 27r) 
and the normalization condition on stationary probability distribution we have the 
expression for stationary current [21] 

A(g) 



J 




-^^''^dq / \{q')e^^'i'Mq' 



G\q) 



dq'dq 



(24) 



From the condition of periodicity it is clear that for the periodic potential and the 
periodic derivative of coupling function with same periodicity V'{q)/[G'^{q)/\{q)] is 
periodic with same periodicity. This makes the effective potential 0(2-7?) equal to 
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Figure 1. Variation of current, J as a function of phase difference, 9 for different 
values Kq and for the parameter set a = 0.5, ksT = 0.1, 7 = 1.0 and D = 1.0. 



zero so that the numerator of Eq. fl2^ reduces to zero. Thus there is no occurrence 
of current for a periodic potential and periodic derivative of couphng function with 
same periodicity and hence there is no violation of second law of thermodynamics. The 
thermodynamic consistency based on symmetry consideration ensures the validity of 
the present formalism. Biittiker [25] have shown that a overdamped particle subjected 
to a drift force field with sinusoidal space dependence and also a sinusoidally modulated 
space dependent diffusion with the same period as the drift experiences a net driving 
force. The resulting current depends on the amplitude of the modulation of diffusion and 
is a periodic function of phase difference between the sinusoidal drift and the sinusoidal 
modulation of the diffusion. 

Let us consider that the particle is moving in a sinusoidal symmetric potential of 
the form 



where Vq is constant and may be taken as barrier height and 6 is the phase factor which 
can be controlled externally. The coupling function is chosen as /(g) = g+a cos g, where 
a is the modulation parameter. We now calculate the current given by Eq. (!24l) . In Fig.l 
the variation of current as a function of phase difference is shown for different values 
of the coupling constant kq. Since is the perturbation parameter in our analysis we 
have kept its maximum value low, i.e., ~ ksT . For the value of the other parameters we 
have chosen a particular set from the complete parameter space. An extensive analysis 
using the full parameter space will be given in our future communication. The current 
shown in Fig.(l) is basically due to the phase difference between the symmetric periodic 
potential and the space dependent diffusion caused by the nonlinear modulation of the 
heat bath by external noise. The current does vanish when the phase difference is either 
zero or integral multiple of vr. When the heat bath is linearly modulated by external noise 
source, it is easy to observe that the effective potential 0(g) is integrable and there will 




(25) 
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Figure 2. Plot of the generalized potential (l){q) as a function of coordinate q for 
different values of Kq and for the parameter set a = 0.5, fcsT = 0.1, 7 = 1.0, D = 1.0 
and 9 = O.Stt. 



be no asymmetry in the effective potential as the noise in the corresponding Langevin 
equation appears additively and the diffusion coefficient becomes space independent. 
Thus when the heat bath is driven nonhnearly by the external noise agency there is a 
net directed motion or phase induced current. This is because of the fact that when 
the external noise drives the heat bath nonlinearly the phase bias gives a tilt to the 
effective potential 0(g) which makes the transition between left to right and right to left 
unequal. In Fig.2 we plot the generalized potential 0(g) for various coupling constant Kq. 
The phase difference (hence the nonlinear driving of the heat bath) breaks the detailed 
balance of the system. When the phase difference is zero or the heat bath is driven 
linearly there is no net drift velocity. Thus, when we drive the heat bath linearly with 
(5-correlated external noise, even in presence of phase difference between V{q) and /'(g) 
there is no net current. For net drift, apart from phase difference nonlinear driving of 
the heat bath is required. This is the central result of this paper. 

In conclusion, we have proposed a new microscopic analysis to study the generation 
of directed motion for a nonlinearly driven heat bath by an external noise. Making 
use of a perturbative treatment we have derived an effective Langevin equation with 
space dependent dissipation and multiplicative noise. Using the corresponding Fokker- 
Planck-Smoluchowski equation with a space dependent diffusion coefficient we have 
checked the thermodynamic consistency condition and have shown that to observe a 
phase induced current one necessarily needs a nonlinear driving of the heat bath by an 
external Gaussian noise. In our future venture in this direction we wish to compare our 
analytical result with stochastic simulation using the complete parameter space. 
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